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ABSTRACT 


The  capabilities  of  the  CRREL  test  basin  to  simulate  dynamic  ice-structure 
interaction  with  scale  model  tests  cover  the  whole  range  of  structures  in 
question.  For  bottom- founded  structure  simulation,  a  test  pile  was  designed 
so  that  its  stiffness,  natural  frequencies  and  modes  and  damping  could  be 
varied.  The  ice  movement  against  the  pile  was  arranged  to  have  constant  accel¬ 
eration  in  order  to  excite  different  modes  with  different  ice  velocities.  The 
flexibility  of  drive  system  caused  jerking  ice  movement  with  low  velocities. 

Analysis  of  the  recorded  ice  forces  and  acceleration  include  the  refinement 
at  measured  ice  forces  by  eliminating  the  response  of  the  measuring  system  itself 
using  dynamic ' equilibrium  or  transfer  function  approach.  The  frequency  contents 
and  the  damping  of  vibrations  are  analyzed  using  Fourier  signal  analyser. 

Scaling  laws  are  discussed  and  it  is  observed  that  full  similitude  cannot 
be  achieved  for  both  ice  interaction  force  and  pile  and  ice  vibrations. 

Results  thus  far  show  similitude  with  in-field  ice  force  histograms:  both 
random  and  saw-tooth  like  repetitious  ice  force  fluctations  appear.  The  change 
from  one  interaction  mode  to  the  other  with  icreasing  drive  velocities  occurs. 
With  high  velocities  relatively  smooth  random  ice  crushing  occurs  and  in  some 
cases  natural  modes  are  stable.  The  greatest  energy  content  of  ice  force  does 
not  always  appear  with  the  natural  frequencies  of  structures. 


1 .  INTRODUCTION 


The  case  of  dynamic  ice  structure  interaction  appears  when  either  the  elas¬ 
tic  deformations  of  the  structure  or  ice,  or  both,  are  affecting  the  ice  failure 
process.  Typically,  this  state  occurs  with  slender  flexible  structures  as 


bottom-founded  piers,  lighthouses  or  moriopods  while  moving  ice  sheet  is  crushing 
against  them.  The  physical  reasoning  lies  on  the  fact  that  the  resistance  of  ice 
is  dependent  on  the  loading  rate  which  on 'its  behalf  depends  on  the  initial  moving 
velocity  of  the  ice  and  on  the  elastic  deformation  rate  of  the  structure  or  the 
ice  at  the  contact  area. 

The  ice  force  fluctuations  and  hence  dynamic  stresses  in  structures  will 
be  considerably  greater  than  in  the  case  where  there  is  no  dynamic  interaction. 
There'  are  two  theories  about  the  origin  of  these  periodically  repeating  force 
fluctuations,  which  in  the  worst  case  may  resonate  with  the  natural  modes  of  the 
structure.  The  first  theory  explains  the  oscillating  ice  force  to  be  a  property 
of  ice  alone,  the  characteristic  of  ice  (Peyton,  1966)  or  the  result  of  ice  to 
break  into  floes  of  certain  size  in  crushing  (Neill  1976).  The  second  theory 
(Maattanen,  1978)  starts  from  the  physical  properties  of  ice-strength  vs.  loading 
rate  and  takes  into  account  also  the  dynamic  behavior  at  the  structure. 

Measurements  for  the  dynamic  ice-structure  interaction  have  been  rather 
limited.  In  most  cases  the  only  objective  has  been  to  measure  the  maximum  ice 
force  values.  Reported  results  this  far  have  indicated  the  existence  of  contin¬ 
uously  repeating  ice  forces  with  frequencies  from  0.5  to  15  Hz  during  in-situ 
measurements.  If  dynamic  interaction  is  not  occuring  the  ice  force  histogram  is 
random  with  smaller  fluctuations.  Only  one  reported  laboratory  test  for  dynamic 
ice  forces  exists  (Peyton,  1966). 

The  drawback  of  field  tests  is  that  the  control  of  the  most  important  para¬ 
meters  ice  strength,  thickness  and  velocity  is  poor  and  chances  to  repeat  the 
test  in  identical  conditions  is  minimal.  In  laboratory  simulations  the  difficul¬ 
ties  are  mainly  in  the  correct  scaling.  However,  the  superior  control  of  main 
parameters,  better  measurement  accuracy  and  chances  for  comparative  tests  more 
than  outbalance  the  scaling  difficulty. 

The  author  of  this  paper  expresses  his  gratitude  to  two  facilities:  to  the 
University  of  Oulu,  Finland,  from  which  he  had  sabbatical  and  to  the  U.S.  Army 
Cold  Regions  Research  and  Engineering  Laboratory  at  which,  during  this  time,  he 
was  able  to  conduct  tests  for  the  dynamic  ice-structure  interaction  in  the  new 
Ice  Engineering  Facility  test  basin.  The  main  objective  of  these  tests  was  to 
measure  the  effects  of  both  the  ice  and  structure  parameters  to  the  rise  of  ice- 

I 

induced  structural  vibrations,  the  ice  force  dependence  on  loading  rate,  amount 
of  ice-induced  positive  and  negative  damping,  and  the  concept  of  dynamically 
stable  structures.  The  following  paper  includes  the  description  of  test  arrange¬ 
ments,  instrumentation,  data  processing  and  of  results  thus  far. 
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'  2-  TEST  SET- VP 

The  CRREL  Ice  Engineering  Facility  teat  basin  is  33.5  ir. 'long,  9.15  m  wide 
and  2.44  in  deep.  The  simulation  of  the  ice  moving  against  a  structure  is  ar¬ 
ranged  so  that  the  test  structure  is  attached  to  the  bottom  of  the  test  basin 
near  the  far  end  and  the  ice  sheet  is  then  nulled  against  it,  Figure  1.  At  the 
rear  end  of  the  basin  there  is  a  slope  and  the  moving  ice  sheet  rises  along  it 
and  falls  down  over  the  edge  to  the  melt  tank. 


Figure  1.  Test  Basin 

The  ice  crystal  size  is  controlled  by  seeding  with  steam  at  15  to  20°C 
whereafter  the  ice  growth  occurred  usually  at  25°C.  Before  tests  the  temper¬ 
ature  is  adjusted  to  predetermined  level.  With  saline  ice  the  temperature  is 
raised  in  such  a  pace  to  get  also  correct  E/o  scaling  (Schwartz,  1975).  The 
walls  of  the  test  basin  are  heated  at  the  water  line  to  keep  ice  from  adfreezing 
to  walls  and  to  guarantee  free  ice  movement  during  the  test.  Open  area  at  the 
walls  is,  depending  on  temperatures  and  the  length  of  warm  up  period,  from  1  to 
10  mm.  Thus,  lateral  confinement  is  minimal  and  the  crushing  situation  is  more 
like  that  of  a  finite  floe. 

After  the  whole  ice  sheet  is  pulled  against  the  structure  and  to  the  melt 
tank,  a  new  freezing  period  can  be  started.  Meanwhile,  the  heat  extracted  from 
making  new  ice  is  used  to  melt  the  previous  ice  in  the  melt  tank.  As  a  curiosity 
it  should  be  mentioned  that  the  whole  heating  of  the  building,  including  offices, 
is  taken  from  the  cooling  of  test  areas.  Capacity  of  the  system  is  to  make  an 
ice  sheet  thicker  than  25  mm  in  a  day.  So  one  test  each  day  can  be  accomplished. 

I 

The  ice  thickness  variation  along  the  length  of  test  basin  has  been  less 
than  +  3%  with  fresh  water  ice  and  less  than  +  10%  with  saline  ice.  The  main 
reason  for  unevenness  is  frost  buildup  on  the  ice  surface.  Ice  upper  surface 
temperature  variations  have  been  usually  less  than  +  0.2®C.  Due  to  difficulties 
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'  in  measuring  ice-sheet  strength  and  modulus  of  elasticity,  most  tests  are  planned 
to  be  comparative:  the  test  basin  length  is  divided  into  two  or  three  parts 
and  a  test  with  only  one  parameter  changed  is  made  subsequently  to  each  of  them. 
As  the  time  required  to  carry  this  through  is  less  than  one  hour,  ice  properties 
can  be  regarded  to  be  the  same  for  each  configuration.  In  later  tests  seismic 
techniques  have  been  utilized  to  measure  the  velocity  of  sound  and  the  modulus 
of  elasticity  of  ice.  The  crushing  strength  is  obtained  in  any  case  as  the  pri¬ 
mary  data. 

The  ice  movement  is  achieved  using  a  7  kW  DC  motor  with  thyristor  drive  to 
turn  reduction  gear  and  winches  which  pull  semi-circular  ice  boom  through  steel 
wire  ropes.  A  semi-circular  ice  boom  was  designed  to  make  it  lightweight  and 
easy  to  handle.  In  addition,  the  circular  boom  exerts  constant  pressure  to  the 
ice  sheet.  The  drive  is  capable  of  ice  forces  up  to  35  kN  and  velocities  from 
0  to  0.20  m/s.  In  a  typical  run  velocity  is  first  increased  from  zero  with 
constant  acceleration  to  maximum  and  then  back  to  zero  similarily  so  that  the 
total  cycle  lasts  about  200  seconds.  This  gives  a  velocity  sweep  to  excite  all 
possible  natural  modes.  As  the  natural  frequencies  are  from  2.5  to  130  Hz,  the 
velocity  is  from  their  point  of  view  pseudo-constant. 

The  major  drawback  of  the  drive  is  the  flexibility  of  steel  wire  ropes. 

This  causes  a  jerking  ice  movement  with  wire  ropes  tightening  and  then  suddenly 
slackening  as  the  ice  fails  against  the  structure.  The  dynamic  interaction  is 
not  between  ice  and  structure  as  intended  but  between  ice  and  drive.  With  higher 
velocities  the  interaction  between  steel  wire  ropes  disappears  but  then  the 
loading  rate  from  the  test  pile  interaction  point  of  view  is  already  far  in  the 
brittle  region.  Installing  hydrodynamic  drag  dampers  in  conjunction  to  the  ice 
boom  stabilizes  a  little  the  jerking  movement  but  does  not.  remove  it  totally  in 
the  low  and  most  interesting  velocity  range.  The  cure  for  this  problem  will  come 
together  with  the  permanent  test  carrier  bridge  which  is  moving  above  the  test 
basin  using  track  and  pinion  traction.  The  test  carrier  can  then  be  used  to  move 
ice  smoothly  on  the  whole  velocity  range. 
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Figure  2.  Test  Pile  and  Supporting  Structure 

The  test  pile,  Figure  2,  is  supported  by  a  rigid  frame  that  is  bolted  to  the 
bottom  of  the  test  basin.  The  dynamic  properties  of  test  piles  can  be  varied 
largely  to  simulate  different  actual  structures.  Its  stiffness  is  varied  by  in¬ 
stalling  the  upper  supporting  clamp  in  different  positions,  its  natural  frequen¬ 
cies  and  modes  can  be  varied  by  varying  the  stiffness  or  by  installing  weights  to 
the  superstructure.  Also  external  damping  elements  can  be  installed.  The  range 
of  spring  constants  for  ice  forces  is  from  0.54  to  8.1  MN/m  and  the  first  natural 
frequency  from  2.0  to  43  Hz  and  the  second  from  7.0  to  80  Hz.  The  waterline  di¬ 
ameters  of  the  pile  have  been  63.5,  101.6  and  177.8  mm.  Also,  different  shapes 
other  than  circular  can  be  easily  installed,  including  inclined  wedges. 

Instrumentation  in  the  test  pile  includes  strain  gages  at  the  support  clamp 
for  the  reaction  of  ice  force  and  at  the  lower  part  of  superstructure  for  the 
elimination  of  superstructure  inertial  effects.  The  test  pile  transducers  were 
calibrated  using  both  static  and  dynamic  known  loadings.  Accelerometers  at  the 
superstructure  are  used  to  measure  the  dynamic  response.  One  accelerometer  is 
installed  on  the  ice  to  measure  ice  jerking. 

Both  analogical  and  digital  recording  is  used.  All  the  signals  are  recorded 
to  instrumentation  tape  recorded  and  simultaneously  strain  gage  signals  and  drive 
RPM  are  digitized  using  a  data-logger.  Digitized  signals  are  stored  on  floppy 
disks.  All  the  instrumentation  is  controlled  by  calculator  during  measurement. 
The  scheme  of  data  recording  and  processing  is  shown  in  Figure  3. 
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Figure  3.  Signal  Recording  and  Processing 
In  a  typical  test  run,  first  the  initial  state  and  a  shunt  calibration  of 
each  data  channel  are  recorded  both  analogically  and  digitally.  Then  the  calcu¬ 
lator  monitors  either  the  drive  RPM  or  the  ice  force  and  whenever  a  predetermined 
level  is  exceeded,  switches  on  the  tape  recorder,  signal  analyzer,  and  starts 
data  logging.  At  the  end  of  the  test  final  state  readings  are  recorded. 

3.  ANALYZING  PROCEDURES 

The  main  interest  in  the  data  processing  is  the  real  ice  force  between  the 
ice  and  the  test  pile.  As  the  direct  force  measurement  is  not  possible  correc¬ 
tions  for  the  dynamic  response  of  the  structure  or  load  cell  has  to  be  calculated. 
It  should  be  noted  that  all  published  dynamic  ice  force  records  thus  far  are  not 
pure  ice  force  records  but  include  also  the  response  effects  of  the  structure 
and/  or  load  cell. 


I  » 


In  digital  signal  processing  the  dynamic  inertia  effects  of  the  superstruc¬ 
ture  are  eliminated  using  the  moment  and  shear  force  values  from  strain  gage 
measurements  just  above  the  waterline.  Although  the  inertial  effects  of  the 
pile  below  waterline  is  not  taken  into  account,  dynamic  calibrations  proved  tV.it 
it  is  possible  to  obtain  the  coefficient  of  correlation  from  0.92  to  0.99  between 
the  calculated  ice  force  and  the  known  excitation  force  with  a  frequency  range 
up  to  65  Hz.  In  the  case  of  resonances  the  coefficient  of  correlation  worsened 
to  0.23  due  to  phase  shifts. 

Outputs  from  the  calculator  are:  corrected  ice  force  histogram,  ice  force 


•vs.  loading  rate  plot,  maximum  and  minimum  values  of  the  ice  forces,  its  mean 
and  standard  deviations. 

In  analogical  signal  processing  in  reality  also  digital  the  correction 
to  ice  force  is  done  using  transfer  function  between  ice  force  and  the  pile 
support  gage.  The  transfer  functions  were  first  measured  using  a  shaker  excita¬ 
tion  with  a  frequency  sweep  from  5  to  150  Hz  but  later  a  simpler  impact  hammer 
method  was  adopted.  In  the  latter  an  impact  through  a  load  cell  is  given  at  the 
point  of  ice  action.  The  impact  gives  an  even  energy  distribution  from  0  to  300 
Hz  and  it  can  be  easily  accomplished  at  the  beginning  of  each  test.  The  signal 
analyzer  then  calculates  the  transfer  function  from  the  signals  of  impactor  load 
cell  and  the  ice  force  reaction  support  gage. 


The  real  ice  force  can  then  be  solved  as  follows.  The  measured  signal  f(t) 
can  be  expanded  into  Fourier  integral  -  extension  of  Fourier  series  expansion  - 
and  in  the  case  of  forced  vibration  it  will  be 

oo 

f(t)  =  -| -  /  H(m)SF(io)e1“tda)  (1) 


where  H(u)  is  the  transfer  function  and  SF(w)  the  linear  spectrum  of  the  forcing 
function.  The  Fourier  transform  for  Equation  1  gives 


H(u)SrM  =  /  f(t)o-ia,tdt  =  Sf(u) 


(2) 


yhich  is  the  linear  spectrum  of  the  measured  signal  f(t).  From  this  Spfui)  can 
be  solved  and  then  the  inverse  Fourier  transform  gives  the  requested  real  ice 
force 


F(t) 


(3) 


Although  the  above  equations  look  complicated  the  signal  analyzer  solves  them 
quickly  using  numerical  fast  Fourier  transform  algorithms.  Same  calculations 
could  be  carried  out  similarly  in  any  computer  for  digitized  data.  As  this 
transfer  function  approach  is  the  only  way  to  obtain  the  real  ice  force,  it  is 
recommended  that  it  is  used  also  in  field  measurements. 


In  addition  to  the  real  ice  force  plots,  the  signal  analyzer  is  used  for 
other  time  and  frequency  domain  analysis:  power  spectras,  cross  spectras,  phase 
information,  and  damping.  The  raw  data  that  has  been  processed  can  then  be 
stored  in  refined  digital  form  to  tape  cassettes  in  the  signal  analyzer. 
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Another  way  for  signal  analyzing  is  comparisons  to  the  predictions  of  mathe¬ 
matical  models.  The  FEM  model  of  the  test  pile  can  be  calibrated  to  give  stiff¬ 
ness  and  the  first  natural  frequency  exactly.  For  dynamic  modc-l  the  dampings 
measured  in  calibrations  can  be  utilized.  Then  by  comparing  the  measured  response 
to  the  predicted  and  the  existence  of  vibrations  to  the  predicted  stability 
criterions,  the  amount  of  ice- induced  positive  or  negative  damping  can  be  figured 
out. 

4.  SCALING  LAWS 

The  scaling  down  of  a  real  structure  for  model  tests  can  be  done  geometri¬ 
cally.  Exact  geometric  scaling  is  not  necessary:  it  is  required  that  the  dy¬ 
namic  behavior  -  naLural  modes  and  frequencies  -  follows  that  of  the  full  scale 
structure.  In  the  case  of  ice-structure  interaction,  it  is  enough  to  scale  only 
those  natural  modes  that  contribute  significantly  to  tiie  deformations  at  the  ice 
action  point.  In  field  observed  range  for  frequencies  is  up  to  15  Hz. 

The  discretized  dynamic  equations  of  motion  of  the  structure  are: 

[k]{6]  +  [d]  {6}  +  [m] { 6 }  =  { F(t) }  (4) 

/ 

where  [k],  [d]  and  [m]  are  the  stiffness,  damping  and  mass  matrices  of  the  struc¬ 
ture,  {6}  nodal  displacement  vector,  16}  and  {6}  velocity  and  acceleration  vector 
respectively,  and  {F(t)}  the  excitation  ice  force.  If  geometric  scaling  factor  n 
and  ice  strength  factor  p  are  applied  and  the  material  of  the  model  and  the 
original  structure  kept  the  same,  then  for  model  with  index  m  the  following 

3 

scaling  terms  are  derived  starting  from  the  beam  bending  stiffness  term  k=  EI/L 
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Similarly,  comparing  term  by  term  in  the  differential  equation,  the  non-diagonal 
terms  for  [k],  [d],  and  [m]  (not  given  in  Equation  5),  car.  be  derived. 

The  in-plane  vibrations  of  the  ice  sheet  are  governed  by  the  differential 
equations: 


— — — — 
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where  u  and  v  are  the  displacements  in  the  x  and  y  directions,  v  the  Poisson's 
ratio,  E  modulus  of  elasticity  and  p  the  density  of  ice.  The  coefficient  on  the 

'i.  2  ' 

right  hand  side  p(l-v  )/ E  is  equal  to  1/c  in  which  c  is  the  plate  compressional 
wave  velocity.  Thus  the  measuring  of  c  by  seismic  tests  is  recommended.  As  can 
be  seen,  the  geometric  and  time  scaling  similarity  as  in  Equation  5  guarantee  the 
similitude  provided  the  ice  properties  are  the  same.  The  ice  grain  size  scaling 
geometrically  is  an  essential  requirement. 


If  ice  strength  is  scaled  by  p  in  Equation  5,  the  displacements  of  the 
structure  are  scaled  np  times.  Then  also  ice  displacements  u  and  v  shall  be  np 
times  from  that  of  full  scale.  The  similitude  then  requires  the  modulus  of 
elasticity  to  remain  the  same.  This  is  usually  impossible,  since  in  practice  if 
ice  strength  is  scaled  down,  the  modulus  of  elasticity  will  change  much  faster. 
This  suggests  that  the  geometric  scaling  n  only  should  be  applied  with  ice  sca¬ 
ling  p  =  1. 


In  dynamic  ice-structure  interaction,  the  appearing  ice  force  depends  on  the 
relative  velocities  between  the  structure  and  the  ice.  While  p  =  1,  this  impor¬ 
tant  parameter  would  be  unchanged  in  the  scaling,  e.g.  vTith  n  =  1:10  and  p  =  1 
the  displacements  of  the  model  and  ice  will  be  one  tenth,  velocities  the  same  and 
accelerations  ten  times  greater  than  with  the  full  scale.  Hov’ever,  the  situation 
is  different  for  the  stress  rate  o.  Tor  ice  moving  against  a  circular  pile  with 
velocity  v  it  is  (Blenkarn,  1970) 


where  o  =  ice  crushing  strength  and  a  =  the  half  of  the  pile  diameter.  To  have 

c  2 
stress  rate  identical  for  both  the  model  and  full  scale  requires  p  /n  =  1.  If  p 

=  1  the  similitude  cannot  be  achieved.  On  the  contrary  if  p  f  1  the  ice  in-plane 

vibration  similitude  cannot  be  obtained.  Thus,  in  either  case  a  compromise  has 

to  be  made.  Thus  far  it  is  unknown  how  the  ice  strength  versus  stress  rate 

changes  from  full  scale  to  model. 


Or.c  in  the  scaling  is  the  floe  size.  In  the  test  basin  only  finite 

size  ice  field  is  possible.  As  the  edges  are  free  due  to  wall  heating,  the 
situation  is  that  of  a  floe  acting  against  the  pile.  Lateral  confinement  does 
not  exist,  but  if  no  longitudinal  cracks  appear,  the  affect  of  boundary  conditions 
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■  will  be  insignif  icant  to  crushing.  The  corresponding  ilex.-  size  in  full  scale 
would  be  typically  greater  than  100  m  by  300  m. 

5.  SOME  RESULTS 

A  typical  ice  force  plot  is  presented  in  app.  1.  At  first  starring  with  low 
velocity,  there  is  interaction  with  the  ice  pulling  steel  wire  ropes  and  the  ice 
force.  After  9  cm/s  this  interaction  is  no  longer  possible  and  then  interaction 
starts  with  the  pile.  In  some  cases,  e.g.  in  app.  2,  while  the  velocity  was 
greater  than  13  cm/s  the  pile  stabilized  and  the  crushing  occurred  with  an  almost 
constant  ice  force.  If  stress  rate  is  calculated  from  Equation  7,  it  appears 

that  in  this  case  minimum  loading  rate  is  with  the  smallest  ice  force  and  highest 

velocity.  This  is  in  contradiction  to  uniaxial  tests,  which  suggests  that  either 

the  Equation  7  is  not  valid  or  that  in  Equation  7  only  the  reference  value  of 

should  be  used  and  not  the  actual  value  for  each  particular  loading  rate. 

The  effect  of  interaction  mode  is  paramount  to  the  ice  forces.  In  all  cases 
the  highest  forces  were  recorded  while  the  interaction  with  the  pulling  system 
was  occurring  in  low  speed.  The  reduction  with  higher  loading  rates  was  consid¬ 
erable,  down  to  40-70%  from  the  maximum  and  as  an  ultimate  to  10%  in  test  27  in 
app.  2.  The  reduction  occurs  both  with  saline  and  fresh  water  ice.  Figure  4 
shows  the  results  of  ice  force  vs.  loading  rate. 
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Figure  4.  Ice  Force  vs.  Loading  Rate 
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Each  individual  point  is  t lie  result  of  3.6  second  period  iron,  the  curve  in  . 

1.  As  can  be  seen  the  maximum  occurs  with  200  rpm*  2.3  cre/s,  and  thereafter  :  •  t: ri 
the  maximum,  average  and  standard  deviation  are  all  decreasing.  The  reduction  in 
the  standard  deviation  is  the  least.  This  type  of  presentation  is  useful  in 
design:  the  probability  of  maximum  ice  force  Lo  become  greater  than  the  average 

plus  three  times  the  standard  deviation  is  very  small.  Thus  an  ice  force  of  (4.6 
+  3.20)kN  =  10.6  kN  extrapolated  to  full  scale  'would  be  an  adequate  design  load. 
As  ice  structure  interaction  is  a  nonlinear  phenomenon  and  as  it  affects  strongly 
the  recorded  ice  force  values,  great  care  should  be  taken  to  observe  interaction 
possibilities. 

Although  the  jerking  ice  movement  prevented  the  ice/test  pile  interaction 
with  velocities  less  than  8  cm/s  the  interaction  with  steel  wire  ropes  gives  some 
information  about  ice  forces  against  ultra  flexible  structures.  The  natural  fre¬ 
quency  of  pulling  systems  varies  from  0.5  to  2.7  Hz  and  if  this  is  changed  to 
full  scale  (1:10)  it  corresponds  to  a  frequency  range  of  0.05  to  0.27  Hz.  The 
cable  tensioning  and  slackening  is  nonlinear  and  it  cannot  sustain  compressive 
loads.  Anyhow,  the  ice  force  repeats  itself  periodically,  app.  I  &  2,  as  long  as 
drive  velocity  is  sufficiently  low.  It  can  be  proved  that  the  vibrating  system 
ice-mass  with  wire  ropes  acting  as  supporting  spring:,  and  excited  with  loading 
rate  dependent  ice  crushing  force  is  dynamically  unstable.  If  the  frequency  of 
ice  force  peaks  is  predicted  using  Equation  8  (Maattanen,  1978) 


max 


(6) 


the  results  will  be  of  correct  order  but  a  little  higher  than  what  it  is  in 
reality.  This  performance  of  Equation  8  coincides  with  another  comparison  of 
Equation  8  to  more  accurate  theoretical  limit  cycle  frequencies.  Due  to  too  many 
and  too  little  known  nonlinear ities  in  the  steel  wire  rope  pulling  system,  no 
mathematical  model  has  been  developed  for  this  case. 


With  increasing  velocity  the  ice  force  excitation  frequency.  Equation  8, 
goes  up  and  after  a  while  interaction  between  steel  wire  ropes  is  no  longer 
possible.  As  this  is  what  happens  with  the  model,  it  supports  the  theory  that 
the  reason  for  ice-induced  vibrations  is  the  result  of  the  ice  strength  depen¬ 
dence  on  the  loading  rate  rather  than  the  result  of  ice  tending  to  break  into 
floes  of  certain  size.  Also,  e.g.  with  300  rptn,  3.4  cm/s,  the  floe  size  should 
be  about  4  cm  instead  of  from  1  to  3  mm  rubble  as  observed  in  crushing.  The 
dynamic  interaction  theory  also  neatly  explains  the  shift  from  the  pulling  system 
interaction  to  the  ice/pile  interaction  with  increasing  velocity,  see  app.  1  &  2, 
and  to  the  stable  state  in  app.  2. 
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Ir.  app.  3  there  arc-  time  expansions  of  the  ice  force  p^ot  in  app.  2.  In  t; 
first  sample  the  drive  velocity  is  low  but  the  energy  stored  in  the  cables  ac¬ 
celerates  the  ice  to  so  high  a  velocity  that  crushing  occurs  with  same  frequenci 
as  with  higher  drive  velocity  in  sample  2.  The  pile  has  natural  frequencies  of 
A3  and  103  Hz.  The  ice  force  excites  the  first  fundamental  mode  slightly  but 
most  of  the  energy  of  ice  force  is  in  the  82  to  87  Hz  range  which  is  lower  than 
the  second  natural  frequency.  With  the  highest  drive  velocity  neither  of  these 
frequencies  is  present,  Lhe  pile  is  stable.  In  1:10  scale  this  situation  would 
simulate  crushing  for  a  structure  having  natural  frequencies  of  A. 3  and  10.3  Hz, 
pile  diameter  1.0  m  and  ice  thickness  59  cm. 


C  S 


The  amount  of  data  and  observations,  and  possibilities  to  simulate  the 
dynamic  ice  structure  interaction  in  laboratory  tests  is  so  great  that  it  is  out 
of  the  scope  of  this  paper.  More  detailed  results  will  be  published  in  a  CRREL 
report  after  the  tests  with  low  speeds  have  been  completed.  In  addition  to 
topics  briefly  treated  here,  it  will  include  data  on:  dynamic  aspect  ratio 
effect,  velocity  dependence  of  crushing  to  buckling  mode  change,  corrections  for 
measured  ice  forces,  existence  of  dynamically  stable  structural  modes,  capa¬ 
bilities  of  predicting  crushing  force  frequencies  and  amplitudes,  ice-induced  po¬ 
sitive  or  negative  damping. 
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